Extraneous Frequencies Generated in Air Carrying Intense 

Sound Waves * 

By A. L. THURAS, R. T. JENKINS and H. T. O'NEIL 

The exact equation of the propagation of plane sound waves in air is not 
linear and consequently harmonics and combination tones are generated. 
The pressure of these extraneous frequencies in terms of the fundamental 
pressure, frequency, and distance from the source has been mathematically 
determined by Rayleigh, Lamb and others. These equations have been 
applied to an exponential horn. 

Measurements of the second harmonic and combination tones have been 
made at various points within a long tube, and in front of an exponential 
horn. Measurements, in general, agree with theory, but the absolute values 
are lower than the calculated values. 

RECENT developments in horn type loud speakers ! for high 
quality reproduction of intense sounds necessitate a consideration 
of the more exact equations of wave motion if distortion due to the 
generation of extraneous frequencies in the air of the horn itself is to be 
avoided. Similar considerations may be of some importance in con- 
nection with the pick-up of intense sounds. 

The propagation of waves of finite displacement has interested 
physicists for more than a century. In 1808 Poisson derived an 
equation which shows that, in general, a sound wave cannot be 
propagated without a change in form and consequent generation of 
additional frequencies. This distortion is caused by the non-linearity 
of air; that is, if equal positive and negative increments of pressure are 
impressed on a mass of air the changes in volume of the mass will not 
be equal; the volume change for the positive pressure will be less 
than the volume change for the equal negative pressure. An idea of 
the nature of the distortion can be obtained from the adiabatic curve 
AB for air as given in the familiar volume pressure indicator diagram 
(Fig. la). The undisturbed pressure and specific volume of air are 
indicated by point P V . Any deviation from the tangent through 
this point causes distortion and consequent generation of extraneous 
frequencies. The theoretical magnitudes of the waves of extraneous 
frequencies are obtained from a solution of the exact differential 
equation of wave propagation in air. The solution shows that the 
pressure of the second harmonic frequency, which is generated in the 
air, increases with the frequency and the magnitude of the fundamental 

* Published in the January 1935 issue of the Jour. Acous. Soc. Am. 
1 E. C. Wente and A. L. Thuras, "Loud Speakers and Microphones," Bell Sys. 
Tech. Jour., 13, 259 (1934). 
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pressure and also with the distance from the sound source. The 
solution also gives the magnitudes of the waves of sum and difference 
frequencies generated when two tones are simultaneously impressed 
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Fig. la — Adiabatic curve for air. 



on the air; these magnitudes also increase with distance from the 
source and with the product of the fundamental pressures and, re- 
spectively, with the sum and difference frequencies. 
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Theory of Propagation of Plane Waves of Finite Amplitude 
The derivation of the exact differential equation for sound wave 
propagation in air involves the continuity equation, Newton's force 
equation and the equation expressing the relation between pressure 
and specific volume in a gas. Since there may be some question as to 
the accurate definition of the density and force in the equation of 
motion a somewhat detailed discussion of this subject will be given. 

Following Rayleigh, 2 let y and y + (dy/dx)dx be the actual distances 
at time t from the plane x = to neighboring layers of air whose un- 
disturbed positions are defined by x and x + dx, respectively, Fig. lb. 
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Fig. lb. 

The displacement corresponding to y is thus £ = y — x and the 
equation of continuity of the fluid is 



p = p (dy/dx)- 1 = Po (l + d£/dx)-\ 



(1) 



where p and p are the densities of the fluid in the disturbed and un- 
disturbed states, respectively. If the effect of viscosity is neglected 
the exact equation of motion of the element of mass p(dy/dx) -dx is 



or 



d 2 y dy, d 2 y dp dy . 

W P dx dX = 6F' Podx= -Ty-fx dx ' 

Po(dH/dt 2 ) = - dp/dx, 



(2) 



p is the pressure at the point y (Fig. lb) which moves with the air 
particle, not the pressure at a fixed point. Except for very large dis- 
placements these pressures are nearly the same. From equations (1) 
and (2) 

PS/dfi = (dp/dp) • (1 + dt/dx)-* ■ (d^/dx 2 ). (3) 

2 Lord Rayleigh, "Theory of Sound," 2nd Ed., Vol. II, p. 31. 
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By virtue of. equation (1), equation (3) is linear in £ only if dp/dp 
= Kp~ 2 or dp/dv = — K, where v = 1/p = specific volume and K is a 
constant. This condition is not satisfied during any ordinary varia- 
tions of state of a gas, but is approximately satisfied when the variations 
are very small. For isothermal changes we have pv = p v and for 
adiabatic changes : 

PlPo = (vo/v) y = (p/po)\ (4) 

where 7 is the ratio of the specific heats and p is the undisturbed 
atmospheric pressure. In either case, for very small variations, the 
pv curve is practically identical with the tangent to the curve, hence 
dp/dv is practically constant (Fig. la). 

From equations (1), (3), (4) we obtain the exact equation of adiabatic 
plane wave motion in a non-viscous fluid : 

32^2 = c \\ + dZ/dxyy-Kd^/dx*), (5) 

where c 2 = ypo/po. This equation is given by Rayleigh. 2 - 8 Rocard * 
was first to call attention to the generation of harmonics in the air 
within an exponential horn. His theoretical solution is based on a 
plane wave equation in which the term d 2 £/d/ 2 was replaced by 

dt 2 dt dx\dt ) K ' 

In support of this substitution Rocard cites Riemann's 6 treatment of 
the problem. However, Riemann's analysis is based on the Eulerian 
form of the hydrodynamical equations whereas equation (5) is derived 
from the Lagrangian equations. (For a comparison of these systems 
of equations see Lamb. 6 ) In the Lagrangian notation d£/d/ and 
d 2 £/67 2 are the exact values of the velocity and acceleration, respectively, 
of the particle whose displacement from its equilibrium position (x) is |. 
It is to be noted that in equation (2) the term p , or undisturbed 
density, does not represent an approximation. 

A rigorous solution of (5) for the displacement £ as an explicit 
function of x and / has not been obtained. As a first approximation to 
equation (5) we take 

dt 2 dx 2 {y + ; dx dx 2 K ' } 

3 Lamb, "Dynamical Theory of Sound," 2nd Ed., p. 182. 

4 Y Rocard, "Sur la propagation des ondes sonores d'amplitude lime, Cotnptes 
rendus, 196, 161 (1933). , 

5 Riemann, " Ueber die Fortpflanzung ebener Luftwellen von endlicner Scnwingung- 
sweitte," Gottingen Abhandlungen, No. 8, 1860. 

6 Lamb, "Hydrodynamics," 6th Ed., Chapter I. 
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This approximation restricts the dilatation d£/dx to values small com- 
pared with unity or the excess pressure to values small compared with 
ypo, but the restriction need not be as severe as would be required for 
the linear approximation: 

d 2 £/dt 2 = c 2 d 2 £/dx 2 . 

By a method of successive approximations, carried to the second 
approximation, Lamb 3 derives the solution of (7) : 

£ = a cos co(t — x/c) -\ ^ ^ a 2 x[l — cos 2co(t — x/c)~\ (8) 

corresponding to a motion £ = a cos cot imposed on the air at x = 0, 
and assuming complete absence of reflection. By virtue of (4) and (1) 
we have for the pressure: 

P= po(l - yd£/dx + ■■■). 

Neglecting the terms of small amplitude, in the region where Awx is 
large compared with the wave-length X we have 



where 



or 



where 



p = Pdc + Pi + Pi, (9) 

Pic = Po - ypo-((y + l)/8) ■(«■/«*)<*», 
pi= — yp -(coa/c) sin co(t — x/c), 
Pi = ypv((y + l)/4)-(wVc»)-a a *sin 2<o(t - x/c), 

pi = 2*Pi cos |>(f - x/c) + tt/2], 

pz = 2*P 2 cos [2«(J - x/c) - tt/2], (10) 

Pi = yp coa/2k, (11) 

Pi ~W 7* T" (12) 

Pi and Pi are thus the r.m.s. fundamental and second harmonic 
pressures, respectively. 

Lamb 3 also gives a solution of (7) for the case when the forced 
motion at x = is: £ = £a cos o>a* + £ B cos co B t. In addition to the 
two fundamentals and two second harmonics, the pressure now in- 
cludes components whose frequencies are, respectively, the sum and 
difference of the two primary frequencies: 

p. = 2>P. cos [(coa + u B )(t - x/c) - tt/2], 
pi = 2*P d cos [(wa - co B )(t - x/c) + tt/2], 
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where 

p _ 7 + 1 PaPr (ma + <»b)x (l3) 

r * 2(2)* ' yp ' c ' 

p _ 7 + 1 PaP b p,i ~ <*b)x , J4 v 

jr " ' 2(2)* 7 po ' " c ' K 

and Pa, Pb are the two r.m.s. fundamental pressures. 

If we extend Lamb's method of solution of equation (7) to the third 
approximation and again consider the case £ = a cos u>t at x = we 
find that in the region kx » 1, the r.m.s. third harmonic pressure is 

In the case of the greatest r.m.s. fundamental pressure used in the 
experiments, Pi = 8000 bars at 600 c.p.s., equation (15) indicates 
that the third harmonic at 400 cm from the source is about 10 db 
below the second harmonic. 

An approximate correction for the effect of attenuation in a tube 
caused by viscosity and heat conduction can be obtained by assuming 
that each of the extraneous frequencies and the fundamental is 
attenuated as if it were the only wave present. Thus the r.m.s. value 
of the fundamental at any point x is assumed to be 

Pi - P e-°i x or dPi/dx = - atPi, 

where P is the r.m.s. value of the fundamental at the point x = 
and «i is the measured attenuation factor for the fundamental. 
If P 2 is the r.m.s. value of the second harmonic at the point x and «2 is 
the measured attenuation factor for the second harmonic in the absence 
of the fundamental, we have by using equation (12): 

dPi/dx = KP\ - aiPi = KP 2 e- 2a i* - a*P a> (16) 

where 

K = y + 1 — - 
A 2(2)* ' yp Q c 

When ai = and ai = 0, equation (16) is equivalent to (12). The 
solution of (16) which is consistent with the fact that the second 
harmonic vanishes at x = is 

P 2 = [KP 2 l(2 ai - a 2 )]|>- a = x - r** m J 
Hence 

P 2 _ 7+ 1 , Po . g R (17) 

Ti 2(2)* y Po C* K ' (17) 
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where 

^ e -(g 2 -a t )x _ e -2a lX ^ ^ 

(2«! - a 2 )* 2 "*" '"' 

For the tube used in these experiments as may be taken as 2 J «i which 
gives 

R = 1 - ai*/2* + •••. 

Similar correction factors were derived for the other extraneous fre- 
quencies measured in the experiments. 

Measurements of Plane Waves of Finite Amplitude 
in a Tube 

The experimental work consisted of measuring the second harmonic 
generated along a tube. Measurements were also made of the sum 
and difference tones when two fundamental frequencies of equal 
pressure were simultaneously impressed on the air of the tube. For 
the fundamental pressures and distances used in the experiments, the 
magnitudes of the other harmonics and higher order sum and difference 
frequencies were probably small. For high fundamental frequencies 
and pressures, however, these other tones are important, since they 
increase more rapidly with frequency and pressure than the second 
harmonic; for instance, the third harmonic pressure increases as the 
square of the fundamental frequency and as the third power of the 
fundamental pressure. 

A sinusoidal displacement, uniform over the cross section, was 
impressed on the air at one end of a long tube. The tube had an inside 
diameter of 3.8 cm and was 1566 cm long. Measurements were made 
in the first 705 cm only and the remainder of the tube was used for 
obtaining a non-reflective termination. 

A search transmitter, comprising a small tube of 0.08 cm inside 
diameter and 7.5 cm long, coupled to a small condenser microphone, 7 
was used for the measurements. Attenuation in this search tube 
was sufficiently high to prevent either overloading the microphone or 
altering the sound wave propagated within the long tube. The search 
transmitter was connected to a stage of amplification so operated as 
to preclude non-linear distortion. This was followed by a band-pass 
filter which selected the frequency desired in the measurements. The 
filter was terminated by a measuring circuit consisting of a high-gain 
amplifier and a vacuum tube voltmeter. A diagram of the arrange- 
ment is shown in Fig. 2. 

7 H. C. Harrison and P. B. Flanders, "An Efficient Miniature Condenser Micro- 
phone System," Bell Sys. Tech. Jour., 11, 451 (1932). 
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To obtain reliable measurements throughout the length of the test 
tube it was necessary to reduce standing waves to a negligible magni- 
tude. This was accomplished by laying a strip of felt in the last 761 cm 
of the tube, terminating the end with an acoustic resistance approxi- 
mately equal to the characteristic impedance of the tube, and carefully 
sealing up all joints along the tube. The pressure variation in the 
standing wave was ± 0.3 db, which corresponds to a reflection coeffi- 
cient of 0.035. 
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Fig. 2 — Apparatus for measuring extraneous frequencies generated in air carrying 
intense sound waves. R, resistance substitute for receiver. 

The oscillator current was supplied to the loud speaker through a 
low-pass filter and the measured harmonic content was found to be 
73 db below the fundamental. Pressure measurements in the tube 
close to the loud speaker indicated that the harmonics generated in the 
measuring circuit and loud speaker were more than 50 db below the 
fundamental pressure at 2000 bars. 

A calibration of the search transmitter was obtained by comparison 
with a small condenser microphone whose diaphragm was exposed 
directly to the sound wave at the same position on the test tube, see 
Fig. 2. The calibrating microphone (C.T., Fig. 2) had been previously 
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calibrated by a thermophone. 8 The measuring circuit following the 
search transmitter was calibrated for each frequency measured by 
introducing the oscillator current into the search transmitter circuit 
through the attenuator (ai, Fig. 2). 

The ratio of the pressure of the frequency generated along the tube 
to the fundamental pressure was measured by the attenuator 0:2, 
Fig. 2, at various holes along the tube in which the search transmitter 
was inserted. If an appreciable fraction of the harmonic generated 
along the tube is reflected at the end of the tube, the magnitude of 
the reflected component near the source may be comparable with or 
larger than the harmonic generated between the source and the point 
in question. This was found to be the case when several measure- 
ments were made near the source over a distance covering a wave- 
length. The measured variation in the total pressure of the harmonic 
over this distance was ± 2.5 db whereas the variation for a funda- 
mental of this frequency, as previously stated, was ± 0.3 db. There- 
fore the measurements close to the receiver may be inaccurate. 

Figure 3 shows the measured pressure ratio of the generated second 
harmonic to the fundamental along the tube and the theoretical curve 
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Fig. 3 — Magnitude of 2nd harmonic vs. distance from source. 

calculated from equation (17). Each of the three experimental points 

plotted at about 45 cm from the receiver is the average pressure ratio 

for a series of readings taken over a distance of a wave-length along 

the tube. 

The measured and theoretical pressure ratios of the second harmonic 

to the fundamental are shown as a function of frequency and pressure 

in Figs. 4 and 5. 

8 L. J. Sivian, "Absolute Calibration of Condenser Transmitter," Bell Sys. Tech. 
Jour., 10, Jan. (1931). 
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Fig. 4 — Variation of 2nd harmonic magnitude with frequency of fundamental. 
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Fig. 5 — Variation of 2nd harmonic pressure with fundamental pressure. 

Figure 6 shows the magnitude of the first order sum and difference 
frequencies along the tube when two frequencies of equal pressure are 
simultaneously impressed on the air in the tube. 
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Exponential Horn Theory 

The second harmonic generated in any short section of a horn is 
approximately the same as that generated in a tube of area equal to 
the mean area of the section of the horn. Therefore, from the tube 
equation and the expression for the change in pressure due to the 
divergence of the horn the magnitude of the generated second harmonic 
pressure at any point along the horn can be obtained. 

The r.m.s. value of a small excess pressure in an exponential horn 
of section 5 = Soe mx is attenuated according to the law 

P = P t e- mi '' 2 or dP/dx = - mP/2, 

where P t is the r.m.s. pressure in the throat of the horn (at x = 0). 
The index of taper m is equal to 4irf c /c where f c is the cut-off frequency 
of the horn. 

From equation (12), the rate at which the r.m.s. value of the second 
harmonic increases along a tube is 



dP 2 7+1 Pi 



a 



dx 2(2)* yp c 



= KP\. 



If it is assumed that the same expression represents the rate of genera- 
tion of second harmonic along a horn, and that both the fundamental 
and second harmonic diverge in the same manner, the complete 
differential equation for Pi becomes 

dp * - KP 2 _ f!t p. = KP* r-m* _ n l Pn 

where Pi and Pn are, respectively, the r.m.s. fundamental pressures 
at the point in question (x) and in the throat of the horn. The solution 
consistent with the condition P* = at x = is 

P 2 = DfiTPf,/(w/2)][>— <' 2 - r— J (18) 

Since Pi = P u e~ mx ' 2 , the ratio of the second harmonic pressure to the 
fundamental pressure at any point x in the horn is thus 

* mXPu .L^2.2gi. £».£*, (19) . 

Pi m/2 2V2 7/> c 

* According to this equation the magnitude of the second harmonic pressure 
generated in the air of a horn is 6 db lower than that given by Rocard's equation 
previously published in the paper "Loud Speakers and Microphones" (Reference 1), 
page 264, equation 1. 
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where 



x ' = (1 - e - m '' i )/(m/2). 



(20) 



Equations (19) and (20) indicate that the second harmonic at the 
mouth of an exponential horn of length x is equivalent to the second 
harmonic at the end of a tube of length x' and of uniform section, 
equal to the area of the throat of the horn. Thus the second harmonic 
in the mouth of a horn having an index of taper m = 0.075 cm" 1 , cut- 
off frequency 200 c.p.s. and length 78 cm is equal to the second har- 
monic at the end of a straight tube 25 cm long and of the same diameter 
as the throat of the horn. 

Exponential Horn Measurements 

Measurements of the output of a horn attached to a moving coil 
receiver were made in an acoustically damped room. The horn had a 
throat diameter of 3.8 cm, a length of 78 cm and a cut-off frequency of 
200 c.p.s. The diaphragm of the loud speaker was coupled to the 
throat of the horn through a straight tube 13 cm long. A filtered 
single frequency tone was impressed on the loud speaker and the sound 
was picked up by a small microphone in front of the horn. The funda- 
mental and second harmonic voltages from the microphone amplifier 
were separated by means of a band-pass filter and measured. 

The approximate acoustic power at the throat of the horn was 
calculated from the known efficiency of the loud speaker and the 
electrical voltage and current supplied. 

The measured and calculated ratios of the second harmonic pressure 
to the fundamental pressure at the mouth of the horn, including the 
effects of generation of second harmonic in both the horn and the 
straight tube coupling the receiver and the throat of the horn, are 
shown in Fig. 7 in terms of the sound output in watts. See equations 
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Fig. 7 — 2nd harmonic generated in an exponential horn VS. sound output. 



EXTRANEOUS FREQUENCIES 



171 



(17) and (19). The attenuation due to viscosity and heat loss in the 
horn has been neglected. 

A number of measurements at various microphone positions in 
front of the horn, shown by the dotted circles, indicate the difficulty 
of obtaining accurate results in a room. The average of the measure- 
ments at a number of random positions in front of the horn at a con- 
stant sound power output and the measurements at a single position 
for various sound outputs gives the plotted curve which is probably 
not greatly different from that which would be obtained in open air. 

Figure 8 shows the measured and calculated ratios of the second 
harmonic pressure to the fundamental pressure for various funda- 
mental frequencies. 
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Fig. 8 — 2nd harmonic generated in an exponential horn vs. fundamental frequency 
(sound output = 10 watts). 



Demonstration of Extraneous Frequencies 

Almost two hundred years ago Sorge, a German organist, and 
Tartini, an Italian violinist, discovered independently, apart from all 
theory, that the union of two loud independent tones produced a 
difference tone. That this is not entirely a subjective tone produced 
by the ear but is actually present in the air was demonstrated by 
others some years later by the use of a tuned resonator. With modern 
apparatus consisting of power amplifiers, oscillators and tuned electro- 
mechanical vibrators it is a relatively simple matter to show not only 
the difference tone but the summation and harmonic tones as well. 

An exponential horn was attached to the open end of the 1566 cm 
tube previously described but with the damping material removed. 
A moving coil microphone placed in front of the horn picked up the 
complex tone produced when two equal pure tones of 600 and 940 
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c.p.s. were impressed on the air. The microphone voltage was 
amplified and impressed on six torsional vibrators tuned to 340, 600, 
940, 1200, 1540 and 1880 c.p.s. Spherical mirrors attached to the 
vibrators produced on a screen bands of light the amplitudes of which 
were approximately proportional to the relative pressures of the 
various frequencies in the complex tone. 

For the higher power inputs to the loud speaker used in the experi- 
ment, the presence of the sum and difference frequencies and the har- 
monic frequencies was easily observed. At these power outputs the 
quality of the sound was very disagreeable and the fundamental tones 
could hardly be distinguished. 

Conclusion 

The theoretical and experimental determinations of the extraneous 
frequency waves generated in the air within a tube are in good agree- 
ment as regards the variation in magnitude with frequency, distance 
from the source and magnitude of the primary tones. The magnitude 
of the second harmonic is very nearly proportional to the distance 
from the source, to the fundamental frequency and to the square of 
the amplitude of the fundamental pressure. When there are two 
primary tones, the extraneous frequencies generated in the air include, 
as well as the harmonics of the primary tones, frequencies which are, 
respectively, the sum and the difference of the primary frequencies 
and also other higher order tones. The magnitudes of the summation 
and difference tones are very nearly proportional to the distance from 
the source, to the product of the magnitudes of the two primary 
pressures, and in each case, to the frequency of the particular com- 
bination tone. As regards the absolute magnitudes of the generated 
tones in the tube all of the measured values are about 3 db lower 
than the theoretical values. 

Good agreement was obtained also between the experimental and 
theoretical determinations of the second harmonic generated in the air 
within an exponential horn, as regards proportionality to the funda- 
mental frequency and power and also as to absolute magnitude. In 
fact the agreement in absolute magnitude was closer for the horn than 
for the tube, but not much significance should be attached to this fact, 
as the horn theory is developed from the theoretical solution for the 
tube and the horn measurements are known to be less reliable than 
the tube measurements. 



